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Abstract 

The purely leptonic decays D s — > rv and D s — > \w are studied in a sample 
of four million hadronic Z decays collected with the ALEPH detector at the 
LEP e + e~ collider from 1991 to 1995. The branching fractions are extracted 
from a combination of two analyses, one optimized to select D s — > rv decays 
with r — > evv or fivv, and the other optimized for D s — > fiv decays. The 
results are used to evaluate the D s decay constant, within the Standard Model: 
/ Ds = [285 ± 19(stat) ± 40(syst)] MeV. 
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1 Introduction 



The leptonic decays of the D s meson are interesting because they open a window onto 
the strong interactions of the constituent quarks of the D s . The absence of strong 
interactions among the final state particles makes the interpretation of the measurements 
particularly straightforward. The information gained is useful in understanding other 
processes involving heavy pseudoscalar mesons. 

The decay D s — > Iv is the second-generation analogue of charged pion decay, 7r + — > Iv. 
The D s — > Iv decay proceeds via the Cabibbo allowed annihilation of the c and s quarks 
in the D s ; the rate depends on the cs wavefunction at zero separation]?] The annihilation 
amplitude is characterized by the D s decay constant /d s , defined as 

ihsPf s) = (D s |c7^75s|0> . (1) 

In the Standard Model the leptonic branching fraction is then given by 

B(B s ^lv) = ^T D Jl\V cs \ 2 m Ds m*(l-^pi , (2) 
8tt V mfyj 

where td s is the mean lifetime of the D s , V ca is the relevant CKM matrix element, m.D s is 
the mass of the D s , and mi is the mass of the lepton. 

It is worthwhile to measure /d s because it characterizes the structure of the D s meson 
and can be calculated in various theoretical models. Lattice QCD is now generally 
considered to be the most reliable method for calculating the pseudoscalar meson decay 
constants; recent results yield a prediction of /d s = 255 ± 30 MeV [Q . An important 
application of the lattice QCD calculations of decay constants is in the evaluation 
of the third-generation Cabibbo-Kobayashi-Maskawa matrix elements. Constraints on 
these elements are obtained from measurements of B°-B° mixing and rely on theoretical 
estimates of /b, which is experimentally inaccessible at present. A measurement of /d s 
therefore serves as a useful check of the theoretical methods. 

The relative branching fractions for D s — > eu, fiv, and tv are known precisely in the 
Standard Model. From Eq. 0, the ratios of branching fractions are -B(D S — > ev)/B(D s — > 
tv) = 2.5 x 10~ 6 and B(D S — ■> fiv)/B(D s — > tv) = 0.103. The helicity suppression in 
these decays leads to an extremely small branching fraction for the electron mode. An 
investigation of the muon and tau modes is presented in this paper. 

The helicity suppression is lifted if a photon is radiated from one of the quarks, giving 
D s — > 7/z/. The analysis of Ref. |2| leads to an estimate of the ratio r = B(D S — > 
7/xi/)/5(D s fiu) = 0.014 to 0.11 with ^(Ds ^eu) = B(D S 7/iz/). The prediction 
in is smaller, r = 0.0059. In either case, the effect of D s — > 7/z/ decays on the present 
investigation would be small or negligible, and no correction is made for these decays. 

The data sample and detector are briefly described in Section |2|. Two analyses are 
then presented, one optimized to select D s — > tv decays (Section |3|) and one optimized to 
select D s — > \iv decays (Section ||). The systematic errors are discussed in Section |5], and 
the combined results and conclusions are given in Section || 

1 Charge-conjugate states are implied throughout this paper. 
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2 Data sample and apparatus 



The analyses are based on a data sample collected with the ALEPH detector at the LEP 
e + e~ collider in 1991 through 1995, at or near the peak of the Z resonance. The sample 
corresponds to 3.97 x 10 6 produced e + e~ — > qq events. 

Detailed descriptions of the ALEPH detector and its performance may be found 
in 0, |J. The tracking system consists of a high-resolution silicon strip vertex detector 
(VDET), a cylindrical drift chamber (the inner tracking chamber or ITC), and a large 
time projection chamber (TPC). The VDET comprises two layers of double-sided silicon 
strip detectors at average radii of 6.3 and 10.8 cm. The spatial resolution for the r-(j) and 
z projections (transverse to and along the beam axis, respectively) is 12 /im at normal 
incidence. The angular coverage is |cos#| < 0.85 for the inner layer and |cos#| < 0.69 for 
the outer layer. The ITC has eight coaxial wire layers at radii from 16 to 26 cm. The 
TPC provides up to 21 three-dimensional coordinates per track at radii between 40 and 
171 cm, as well as measurements of the specific ionization energy loss (dE/dx) of charged 
particles. The tracking detectors are contained within a superconducting solenoid, which 
produces an axial magnetic field of 1.5 T. Charged tracks measured in the VDET-ITC- 
TPC system are reconstructed with a momentum resolution of Ap/p = 6 x 10~ 4 p t ©0.005 
(p t in GeV/c). An impact parameter resolution of 22 /xm in the r-0 plane is achieved for 
muons from Z — ► having at least one VDET r-<p hit. These performance figures 

reflect the improvements obtained from the 1998 reprocessing of the ALEPH LEP1 data 
set with enhanced event reconstruction algorithms [/j. 

Surrounding the TPC is an electromagnetic calorimeter (ECAL), a lead/ wire- chamber 
sandwich operated in proportional mode. The calorimeter is read out via projective towers 
subtending typically 0.9° x 0.9° which sum the deposited energy in three sections in depth. 
Beyond the ECAL lies the solenoid, followed by a hadron calorimeter (HCAL), which uses 
the iron return yoke as absorber and has an average depth of 1.50 m. Hadronic showers 
are sampled by 23 planes of streamer tubes, providing a digital hit pattern and inducing 
an analog signal on pads arranged in projective towers. The HCAL is used in combination 
with two layers of muon chambers outside the magnet for /x identification. 

The measurements of charged particle momenta and of energy depositions in the 
calorimeters, combined with the identification of photons, electrons, and muons, are used 
to produce a list of charged and neutral energy flow particles in each event |J. 

3 D s — ► tv analysis 

The first analysis || is designed to search for the decay chain e + e~ — ► cc, c — ► D s — > tv, 
with t — > qvv or r — > {ivv. The electron and muon channels are treated separately. As 
the final state under study contains three neutrinos, the signature of these decays is an 
identified lepton and large missing energy in one hemisphere of the event. 

Although the e + e~ — > bb event sample contains more D s mesons than the cc sample, 
the D s — > Iv signal is more difficult to isolate in bb events due to the softer D s spectrum 
and the large number of semileptonic b decays. Cuts are therefore applied to remove bb 
events, and the analyses are optimized to select D s decays in cc events. 
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3.1 Event selection 



Hadronic Z decays are preselected using charged tracks || . Backgrounds from two-photon 
interactions and dilepton events are reduced to a negligible level by means of additional 
cuts on the numbers of reconstructed charged and neutral particles. The event thrust axis 
is required to satisfy | cos thrust I < 0.8 to select events within the acceptance of the vertex 
detector. Each event is then separated into two hemispheres by the plane perpendicular 
to the thrust axis. The total energy of each hemisphere is calculated and the hemisphere 
with larger missing energy is selected if it contains an identified lepton (e or //). Electron 
identification is based on the shower shape in the electromagnetic calorimeter and the 
ionization in the time projection chamber; muon identification makes use of the digital 



pattern information in the hadron calorimeter and hits in the muon chambers [Kj. If 
more than one lepton is present in the selected hemisphere, the one with the highest 
momentum is considered. The missing energy in the hemisphere is required to be greater 
than 5 GeV to reduce the background. The hemisphere invariant masses are taken into 
account in the calculation of the expected hemisphere energies [|ll[ . 



To further reduce the background from bb events, cuts are applied to a set of 
hemisphere variables based on the pseudorapidities and impact parameters of the charged 
tracks. The energy flow particles are clustered into jets by means of the JADE algorithm 
with y cnt = 0.003. The pseudorapidity r\ of a particle is then defined as 77 = — lntan(a/2), 
where a is the laboratory angle between the particle and the nearest jet axis. The 
reconstructed particles in cc events tend to have larger pseudorapidities than those in 
bb events because of the lower mass of c hadrons compared to b hadrons. An existing 
ALEPH lifetime-based b tag algorithm |12| was modified to include a dependence on the 
pseudorapidity of the charged tracks. The tracks in each hemisphere are divided into two 
groups, one with 77 > 5.1 and one with 77 < 5.1. For each group, the confidence level for 
the hypothesis that all tracks originate from the Z production point is then calculated. 
The selected lepton track is not included in this calculation. To reject bb events, a cut is 
made on the confidence level for the low-77 group in each hemisphere. A further rejection 
of cc background events^ is obtained by cutting on the confidence level for the high-// 
group in the hemisphere containing the lepton; with the lepton excluded, signal events 
have only fragmentation tracks in the lepton hemisphere and can be distinguished from 
cc background events. 

The momentum and energy of the D s candidate are reconstructed from the observed 
charged and neutral particles in the event by excluding the lepton I and applying four- 
momentum conservation: 

This procedure is based on the assumption that no neutrinos are produced in the 
hemisphere opposite the D s candidate. A kinematic fit is performed in which the energies 
of all reconstructed particles (except the lepton) are varied such that the constraint 

2 Throughout this paper, light quark-antiquark (denoted "uds"), cc, and bb background events are 
defined to be those events containing no D s — * Iv or D + — > Iv decays. 
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[E^ — P^ c 2 ] 1 / 2 = M Ds c 2 = 1.968 GeV is satisfied. The energy resolutions for charged, 
neutral electromagnetic, and neutral hadronic particles are parametrized from simulated 
events. The kinematic fit improves the energy resolution of D s candidates from 6.7 GeV 
to 3.5 GeV. The background is further reduced by the requirement that the fitted D s 
energy be greater than 22.5 GeV. The selection efficiency for cc — > D s — > tv events is 
2.5% (3.3%) in the electron (muon) channel [including the factor of B(r — > Ivv)}. This 
procedure selects 3956 and 6637 events in the electron and muon channels, respectively. 



3.2 Linear discriminant analysis 

A linear discriminant analysis is performed on Monte Carlo events to search for an optimal 
linear combination of event variables so that maximum discrimination between signal 
and background is achieved. Two discriminating variables, U c and Ub, are created to 
distinguish the cc — > D s — > tv signal events from cc and bb backgrounds, respectively; 
the backgrounds arise mainly from semileptonic c and b decays. Nine variables are selected 
to form U c and ten for Ub- The variables with the greatest discrimination power include 
the fitted D s momentum, the angle between the lepton and the D s direction in the D s rest 
frame, several b tag variables, and the lepton transverse momentum with respect to the 
nearest jet axis. The definitions of U c and Ub are optimized independently in the electron 
and muon channels. The distributions of Ub versus U c are shown in Figs. |T] and ^. 

The branching fraction B(D S — > tv) is extracted from a maximum likelihood fit to the 
unbinned Ub versus U c distribution in each channel. The two-dimensional fitting function 
consists of one signal and three background (uds, cc, and bb) components. Although the 
analysis is optimized to select D s — > tv decays, some D s — * \w decays are also accepted. 
Because these two decay rates are both proportional to , both types of decays are 
considered to be part of the signal in the fit to the data. Leptonic D + decays also pass 
the selection and are included in the signal. The production of D s and D + mesons in the 
remaining bb events must be taken into account as well. The signal component of the 
fitting function is therefore taken to be the combination of D s and D + decays to tv and 
jiv in cc and bb events. The relative normalizations of these eight contributions are fixed 



according to Eq. ^| and the measured charm production rates [|T3j, [TJ[]; /d s //d = 1.11 



+0.06 
-0.05 



from lattice QCD calculations fL5| ; R c and from the Standard Model fit in ||13|| ; \V C 



cs I 



and \V c< i\ from fl6fl ; the D s and D + lifetimes from [ lql ; and particle masses from [[U| 



The normalizations are shown in Table |I[ Leptonic D + decays are Cabibbo suppressed 
and contribute less than 8% of the total signal. Each component of the fitting function 
is parametrized as the sum of up to five two-dimensional correlated Gaussian functions 
from a fit to simulated events. 



A program based on JETSET 7.4 |I7j and tuned to ALEPH data [Fj| was used to 
generate the Monte Carlo sample of 15 million e + e~ — > qq background events (including 
3.7 million cc and 6.6 million bb events) and 370 000 events containing D s and D + 
leptonic decays. The polarization of the r is taken into account in the simulation of 
(pseudoscalar) — > tv decays. 

The free parameters in the fit to the data are the numbers of events in the signal 
and the three background components. The fitting procedure is tested with Monte Carlo 
events to be free of bias. 
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(a) Signal: cc -> D s — > TV 
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(c) Background: Z — > cc 
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(b) Background: uds events 
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(d) Background: Z — > bb 
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Figure 1: Monte Carlo Ub versus U c distributions in the electron channel of the rv 
analysis, after all cuts: (a) signal, cc — > D s — > ru; (b) uds background; (c) cc background; 
(d) bb background. The area of the square in each bin is proportional to the number of 
entries. The distributions are shown with arbitrary normalizations. 
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ALEPH 

(a) Signal: cc — > D — > TV 
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(c) Background: Z — > cc 
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(b) Background: uds events 
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(d) Background: Z -> bb 
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Figure 2: Monte Carlo Ub versus U c distributions in the muon channel of the rv analysis, 
after all cuts: (a) signal, cc — > D s — > tv\ (b) uds background; (c) cc background; (d) bb 
background. The distributions are shown with arbitrary normalizations. 
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Table 1: Relative normalizations of the contributions to the signal in the D s 



tv analysis. 



Fraction of signal (%) 
Source e channel /i channel 



cc -»• D s -> tv 76.9 52.4 

cc->D s ->//z/ 0.5 26.4 

bb -> D s -> rz/ 17.4 10.7 

bb -> D s -> /iz/ 1.2 3.1 

cc -> D+ -> ri/ 3.4 2.3 

cc -> D+ -> /xi/ 0.1 4.4 

bb -> D+ -> n/ 0.4 0.4 

bb -» D+ -» ^ 0.2 04 

Total 100.0 100.0 



Table 2: Fitted numbers of signal and background events in data in the D s — > tv analysis. 
Only the statistical uncertainties from the fit are shown. 



Component 


e channel 


/i channel 


Signal 


306 ± 62 


575 ± 84 


uds background 


111 ± 56 


455 ± 139 


cc background 


2310 ±101 


3750 ± 182 


bb background 


1228 ± 56 


1857 ± 74 



3.3 Results 

The projections of the two-dimensional fits to the data are shown in Figs. |] and f|. The 
fitted numbers of signal and background events are listed in Table ^. The fitted branching 
fractions from the electron and muon channels are -B(D S — ► tv) = [5.86 ± 1.18(stat)]% 
and [5.78 ± 0.85(stat)]%, respectively. 

The goodness of fit is evaluated by comparing the likelihood value obtained from 
each fit to the data with the distribution of likelihoods obtained from fits to many toy 
Monte Carlo samples. The fitting functions are used to define the parent distributions 
in generating these samples. The fit confidence levels are calculated to be 83% in the 
electron channel and 81% in the muon channel; the Ub versus U c distributions in the 
data are well described by the fitting functions. The same projections are also shown in 
Figs. H and f| after subtraction of the fitted background components. The distribution of 
the excess events is consistent in shape with the Monte Carlo prediction for the signal. 

It is desirable to visualize the results of the two-dimensional fits in one-dimensional 
distributions. For this purpose, the purity P of each event is calculated. The purity of 
event i is defined as 

where and Tjj^ are the values of the discriminant variables for event i, and S and 
B represent the distributions of the signal and background, normalized according to the 



7 




Figure 3: Projections of the fit to the data in the electron channel of the D s — > tv 
analysis: (a) U c distribution; (b) Ub distribution. The data are shown by the squares 
with error bars. The solid curves are the fitted distributions, while the dashed curves 
show the signal contributions. The same distributions are shown in (c) and (d), after 
subtraction of the fitted backgrounds. 
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Figure 4: Projections of the fit to the data in the muon channel of the D s — > tv analysis: 
(a) U c distribution; (b) Ub distribution. The data are shown by the squares with error 
bars. The solid curves are the fitted distributions, while the dashed curves show the signal 
contributions. The same distributions are shown in (c) and (d), after subtraction of the 
fitted backgrounds. 
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results of the fit. The distributions of P are quite different for signal and background 
events. The purity distributions from the electron and muon channels are shown in Figs. |5]a 
and b; the background-subtracted distributions are given in Figs. |5|c and d. These figures 
show that the shapes of the observed purity distributions agree with the expectations. In 
particular, the shape of the data distribution in each channel after background subtraction 
is consistent with the simulated signal and inconsistent with the simulated background. 



4 D s — > fits analysis 



The second analysis |19j is optimized to select e + e 
D s -> pus. 



cc events containing the decay 



4.1 Event selection 

The same preselection of hadronic Z decays is performed as in the rv analysis. A cut 
on the number of reconstructed charged tracks coming from the area of the interaction 
point is applied to reduce the background from dilepton events. The event thrust axis is 



required to satisfy |cos 6Wust | < 0.8. A loose muon identification algorithm |19|, requiring 
either muon chamber hits or a muon- like digital pattern in the HCAL, is used to select 
muon candidates. If more than one muon is present in the event, the one with highest 
momentum is selected. 

A kinematic fit is then performed in order to improve the resolution on the missing 
momentum, assumed to arise from the undetected neutrino from D s — > [iv. In this fit, 
the missing mass is constrained to be zero and the D s displacement direction (from the 
event primary vertex to the D s decay vertex, i.e., an unknown point on the muon track) is 
constrained to be parallel to the D s momentum direction. This constraint is equivalent to 
requiring that the reconstructed neutrino direction be parallel to a plane containing the 
primary vertex and the muon track. The energies of the reconstructed charged and neutral 
particles are varied in the fit, while their directions are held constant. The muon track 
impact parameters and the coordinates of the event primary vertex are also allowed to 
vary within their uncertainties. The kinematic fit improves the resolution on the neutrino 
energy from 5.7 GeV to 2.8 GeV and on the neutrino direction from 8.9° to 6.2°. 

The bb background and some of the cc background are further reduced by cuts on 
the track pseudorapidities and impact parameters, as described in Section |3TT| . Finally, 
a hard cut is made on the fitted energy of the D s candidate, + E u > 28 GeV. This 
procedure selects 7164 events in the data. 



4.2 Linear discriminant analysis 

Linear discriminant variables U c and Ub are optimized to distinguish cc — > D s — > fiu signal 
events from cc and bb backgrounds, respectively. Each of the linear discriminant variables 
is a combination of seven event variables. The variables with the greatest discrimination 
power include the muon momentum, the fitted neutrino momentum, and several b- and 
c-tag neural network outputs. The Ub versus U c distributions for signal and background 
events are shown in Fig. || 
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Figure 5: Distributions of purity P, in the (a) electron channel and (b) muon channel 
of the D s — > tv analysis. The data are shown by the squares with error bars. The 
histograms show the Monte Carlo distributions for the background (dotted) and for the 
signal plus background (solid). The Monte Carlo histograms are normalized according to 
the results of the fits to the data. The plots in (c) and (d) show the same distributions 
after subtraction of the fitted backgrounds; the squares are the data and the histograms 
are the Monte Carlo signal distributions. The error bars in all four plots reflect the 
statistical errors in the data. 
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(a) Signal: cc — > D s — > uv 



-1 



□ 




□ 


□ □ ■ 


. = □ □ 


□□□□ □ n = ■ 


r □ □□ 


□□□□□ □ □ = 


■ ■ □ □□ 


Mjnn □ n 


■ □ □□ 


□□□□□ □ - 


■ ° □ □ 


□□□□□ □ □ ■ 


■ □ □ □ 


□□□ □ □ n □ - 


- □ □ □:□ □ □ □ □ - ■ 


= □ □ 
□ □ 




□;□□□□» ■ 

□ : □ □ □ □ 

□ : □ □ 

i 


-1 


1 



(c) Background: Z — » cc 
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(b) Signal: cc -> D -> TV 
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(d) Background: Z — » bb 
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Figure 6: Monte Carlo Ub versus U c distributions in the D s — > [w analysis, after all 
cuts: (a) signal, cc — > D s — > fiu; (b) cc — > D s — > tv, which is also included in the 
signal component of the fitting function; (c) cc background; (d) bb background. The 
distributions are shown with arbitrary normalizations. 
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Table 3: Relative normalizations of the contributions to the signal in the D s — > fiv analysis. 



Source Fraction (%) 

cc - 
cc - 
bb 
bb 
cc - 
cc - 

bb 
bb 

Total 



n 


TV 


^ty .0 


D s - 


-> (IV 


26.7 


D s - 


-> TV 


12.6 


D s - 


-> /XI/ 


3.5 


D+ - 


-> TV 


2.2 


D+ - 


-> /XV 


4.3 


D+ 


— > rv 


0.6 


D+ 


— > /xz/ 


0.5 


100.0 



4.3 Results 

The numbers of signal and background events are extracted from the data by means of a 
maximum likelihood fit to the (U c , Ub, M^-J) distribution, where is the invariant mass 
of the D s candidate after the kinematic fit. The three-dimensional region —0.6 < U c < 1.0, 
-0.8 < U b < 1.4, and < 5 GeV/c 2 is divided into 6x6x11 bins in the fit. The fitting 
functions are constructed from simulated events and a procedure is applied to smooth the 
distributions to reduce the bias that results from limited Monte Carlo statistics near the 
edges of the space. 

As in the D s — > tv analysis, the signal component consists of eight contributions from 
D s and D + decays to tv and \w in cc and bb events; the relative normalizations of these 
contributions are again fixed according to the Standard Model expectations (Table |3|). 
Although the \xv invariant mass provides some separation between D s — > \xv decays and 
the backgrounds, the D s — > fiv decay mode comprises only 30% of the signal and cannot 
be clearly isolated in a fit to the (U c , Ub, M^ v ) distribution. With the Standard Model 
constraint -B(D S — > fiv)/B(D s — > tv) = 0.103 imposed in the fit, this analysis nevertheless 
yields additional information on /d s . (Fits in which this constraint is not imposed are 
described in Section |y|) 

The fitted numbers of events are 553±93 signal events, 166±47 uds background events, 
1251 ± 71 cc background events, and 1291 ± 62 bb background events. Figure [7| shows 
the fit projection on the axis. After correction for a bias of +8.2% in the smoothing 
and fitting procedure, the fit result corresponds to -B(D S — ► \iv) = [0.68 ± 0.11(stat)]%. 
The goodness of fit is characterized by a confidence level of 69%. 

The purity distribution from the D s — > fiv analysis is shown in Fig. |8]. The purity is 
defined as in Section |3~3| , except in this case the signal and background distributions S 
and B are binned functions of U c , Ub, and M^. Again the signal and background events 
have distinct distributions, and the shape of the background-subtracted data distribution 
agrees with that of the simulated signal. 



5 Systematic errors 
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Figure 7: (a) M^ v projection of the fit to the data in the D s — » \iv analysis. The data 
are shown by the squares with error bars. The solid curve is the fitted distribution. The 
dashed curve shows the signal contribution, while the dotted curves indicate the three 
background contributions. The data and fitted signal distributions are shown in (b), after 
subtraction of the fitted backgrounds. 
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Figure 8: (a) Distribution of purity P in the D s — > jj,u analysis. The data are shown 
by the squares with error bars. The histograms show the Monte Carlo distributions for 
the background (dotted) and for the signal plus background (solid). The Monte Carlo 
histograms are normalized according to the results of the fit to the data, (b) Purity 
distribution after subtraction of the fitted backgrounds; the squares are the data and the 
histogram is the Monte Carlo signal distribution. The errors in both plots reflect the 
statistical errors in the data. 
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Table 4: Relative systematic uncertainties on B(D S — > Zz/). 



Uncertainty (%) 



Source 


D s -> ri/ (e) 


D s -> ti/ (/^) 


D s -> /xz/ 


Charm hadron production 


27.0 


22.8 


19.6 


c fragmentation 


16.7 


14.8 


12.1 


b fragmentation 


9.0 


5.9 


4.1 


Lepton spectra in b and c decays 


5.5 


3.6 


3.1 


Lepton rates in bb events 


0.3 


0.1 


1.7 


Detector resolution 


13.6 


12.7 


4.0 


Monte Carlo statistics 


7.2 


6.0 


8.8 


Bias in fitting procedure 






4.1 


D + — ► K°yU + z/ form factor 


0.5 


0.5 


1.5 


/d s / /d 


0.2 


0.7 


1.4 


Total 


36.8 


31.4 


26.0 



5.1 Evaluation of systematic uncertainties 

In both analyses, the shapes of the signal and background fitting functions are taken from 
simulated events. The physics parameters that can affect these shapes are taken into 
account in the systematic uncertainties on the branching fractions. 

The sources of systematic error are summarized in Table f|. The magnitudes of the 
systematic errors are estimated by measuring the branching fractions from a Monte Carlo 
sample, reweighting the Monte Carlo events to simulate a variation in a parameter, then 
measuring the branching fractions again. In each case the quoted systematic uncertainty 
is taken to be the quadratic sum of the observed shift and the statistical uncertainty on 
the shift. 

The uncertainties on the production rates of D + , D°, D s , and charm baryons 



from []13| , |14| are taken into account, including correlations. The D s production rates have 
large uncertainties due mainly to the large uncertainty on the absolute scale of branching 
fractions for hadronic D s decays; these parameters give the largest contributions to the 
systematic error on the leptonic branching fractions because they directly govern the 
number of produced D s mesons in the data sample. 

The Monte Carlo events are weighted to match the mean scaled energies (xe{c)) = 
0.484 ± 0.008 and (x E (h)} = 0.702 ± 0.008 of charm and b hadrons given in 0; the 
uncertainties on those fragmentation parameters are propagated to the leptonic branching 
fractions. 

The shapes of the lepton energy spectra in b and c decays are varied according to the 
prescription given in [^TJ. The b — > I and b — ► c — > I fractions and uncertainties are taken 
from |13| . 



The detector resolution on missing energy and momentum is dominated by the errors 
on neutral particle energies. The energy resolutions for neutral energy flow particles are 
studied in samples of hadronic Z decays in data and Monte Carlo. The electromagnetic and 
hadronic energy resolution parameters for each sample are estimated by means of a global 



16 



fit in which the total energy in each event is constrained to y/s and the total momentum is 
constrained to zero. Although the four-momentum carried away by neutrinos and objects 
at low angles is not taken into account in this procedure, a useful comparison of the energy 
resolutions in data and Monte Carlo can nevertheless be made. A discrepancy of less than 
4% in the neutral particle energy resolution is observed, and the effect on the branching 
fractions is estimated by further smearing the neutral energies in the Monte Carlo events 
that are used to build the fitting functions. 

The Monte Carlo statistical uncertainty is estimated by generating many toy Monte 
Carlo samples, using the fitting functions obtained from the full simulation as the parent 
distributions. This procedure reveals a +8.2% bias in the smoothing and fitting procedure 
in the D s — > \xv analysis, and a systematic uncertainty equal to half the bias is assigned. As 
a cross check, a fit to the data in this channel is also performed with an algorithm in which 



the statistical fluctuations in the Monte Carlo distributions are taken into account | 2~2| ; the 
fitted number of signal events is within 2% of that obtained with the standard program, 
after bias corrections. 

Other small effects considered are the form factor in D + — > K°fi + v decays |23[ (an 
important source of background) and the ratio of decay constants fnj Jd from (TBI . 



5.2 Cross check with D s — > (pn decays 

A study of D s — > <pn decays is performed as a cross check of the signal efficiencies in the 
D s — > Iv analyses. Candidate D s — > <f)n decays with <p — > K + K~ are first selected from the 
full e + e~ — > qq sample. The dE/dx measurements in the TPC are used to discriminate 
between pions and kaons. Cuts are also applied on the kaon and pion momenta, the K + K~ 
invariant mass, and the decay angle of the <fi. The same selection is applied to simulated 
qq events containing D s — > <f)ir decays. 

The KK7T candidates are divided into seven bins of xe = Ekk-k / -E-beam, from 0.3 to 
1. The number of D s mesons in each xe bin in the data is evaluated by means of a fit 
to the reconstructed KKn invariant mass distribution for the candidates in that bin. In 
this fit the signal is described as the sum of two Gaussian functions, and a second pair 
of Gaussians is included for the D + — > 07r contribution; a polynomial function is used for 
the background. 

In order to measure the efficiencies of the D s — > Iv selections, the pion in each KKn 
combination is treated as the lepton candidate; the kaons are removed from the event and 
correspond to the missing neutrino(s). After the selection cuts are made, the number of 
D s mesons in each xe bin is again extracted from the KKn invariant mass distribution. 
The selection efficiency in each bin is then evaluated in data and Monte Carlo. Some 
cuts, most notably the lepton identification cuts, cannot be studied with this method. 
The D s — > tv and D s — > fiv analyses are checked separately. 

The following aspects of the D s — > tv analysis are included in the 07r cross check: 
the cut on missing energy, the bb and dilepton rejection cuts, the kinematic fit, and the 
cut on the fitted D s energy. With the kaon candidates excluded from the event, these 
operations can be accurately duplicated in the 07r events. The overall efficiency of these 
cuts as a function of xe, in data and Monte Carlo, is plotted in Fig. |]. The signal 
efficiency (apart from lepton identification) is verified to be correctly simulated within 
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Figure 9: Efficiency versus Xe in data (squares with error bars) and Monte Carlo 
(histogram) for the D s — > tu selection cuts listed in the text, as measured from D s — > (fin 
decays. The cut on the fitted D s energy, _E Ds > 22.5 GeV, corresponds to x s >0.49. 
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Figure 10: Background-corrected distribution of the D s mass (calculated from the n and 
the missing K + K~) after the kinematic fit, in the <f)ir cross check of the D s — > fiu analysis. 
The data are shown by the squares with error bars, while the histogram shows the Monte 
Carlo distribution. The number of D s mesons in each mass bin in the data is obtained 
from a fit to the corresponding distribution of the fully reconstructed KK7T invariant mass. 
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Table 5: Fitted background events relative to the expected numbers in the two channels 
of the D s — > tv analysis and in the D s — ► fiv analysis. The first uncertainty given in each 
case is statistical and the second systematic. 





Fitted / Expected Background 


Background 


D s -> tv (e) D s -> tv (//) D s -> \w 


uds 


0.81 ±0.40 ±0.58 1.00 ±0.31 ±0.42 0.84 ± 0.22 ± 0.17 


cc 


0.92 ±0.04 ±0.06 0.93 ±0.04 ±0.07 0.96 ± 0.05 ± 0.08 


bb 


0.92 ±0.05 ±0.05 1.06 ±0.04 ±0.05 0.98 ± 0.05 ± 0.06 



the statistical uncertainty of the test (7% relative uncertainty). 

A similar 07r cross check is performed for the D s — > \iv analysis. The main differences 
with respect to the D s — > tv case are (1) the fitted primary vertex and the impact 
parameters of the pion (serving as the muon) are now involved in the kinematic fit, and 
(2) it is possible to check the resolution on the reconstructed D s mass. For purposes of 
this test, the missing mass is constrained in the kinematic fit to equal the reconstructed 
K + K~ invariant mass on an event-by-event basis. The measured efficiencies in data and 
Monte Carlo are again found to be in agreement. The distributions of reconstructed D s 
mass in data and Monte Carlo are plotted in Fig. [LC]; the simulated resolution is found to 
be consistent with that observed in the data. 



5.3 Additional checks 

Additional checks, described in this section, are performed in order to further investigate 
the accuracy of the Monte Carlo simulation. 

The fitted numbers of uds, cc, and bb background events are compared with the Monte 
Carlo predictions in Table [|. The normalization of the predicted numbers is based on 
the total number of produced e + e~ — > qq events in the data sample. The predictions for 
D s — ► fiv are corrected for biases in the fitting procedure. The data and Monte Carlo 
are in agreement. A \ 2 value is computed to characterize this agreement in the D s — ► fiv 
analysis, taking into account the correlations in the statistical and systematic errors on 
the normalizations of the three backgrounds. The result is \ 2 = 0-74 for three degrees of 
freedom (CL = 87%). 

The distributions of the event variables used to construct the three sets of linear 
discriminant variables U c and Ub are compared in data and Monte Carlo to ensure that 
the simulation is reliable. The Monte Carlo plots are normalized according to the results 
of the fits to the data. No significant discrepancies are found. 

Another cross check on the accuracy of the simulation is made by comparing the 
distributions of certain quantities in data and Monte Carlo in different regions of the 
(U c , Ub) space. In each analysis, four disjoint regions in that space are defined with 
different relative abundances of the signal and the three backgrounds. The Monte Carlo 
normalization is again based on the results of the fits to the data. The quantities studied 
in the D s — > tv analysis include the fitted D s momentum, the angle between the lepton 
and the D s momentum in the D s rest frame, and b tag variables; in the D s — > fiv analysis 
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the missing energy before the kinematic fit, the fiv invariant mass after the kinematic 
fit, and the fitted proper D s decay length were considered. The simulation is found to 
be consistent with the data in all cases and in all regions of the (U c , Ub) space. This 
check demonstrates that the properties of the backgrounds are correctly reproduced in 
the Monte Carlo. 

The fit to the data in the D s — > jj,u analysis was repeated without the Standard Model 
constraint on B(D S — > fiv)/B{p a — > tv). For this calculation, the signal contributions 
from D s and D + decays to fiv were allowed to vary separately from the tv contributions. 
The fitted contribution from \xv relative to tv is found to be somewhat smaller than 
expected; Monte Carlo experiments show that the unconstrained fit is expected to give a 
Iw/tv ratio smaller than the observed value 5% of the time when the events are generated 
according to the Standard Model ratio. The fitted [iv/tv ratio is not sensitive to the 
assumed charm production rates or fragmentation parameters. 

Increased sensitivity to B(D S — > jiv) is achieved if the D s — > \w analysis is coupled to 
the statistically independent electron channel of the D s — > tv analysis. The fit described 
in the preceding paragraph is repeated with a modified likelihood function in which 
B(D S — > tv) is constrained according to the value and uncertainty obtained from the 
electron channel; after correction for the calculated —1% and +14% biases, respectively, 
the results are 

B(D a ^fiu) = (0.47 ±0.25)% 
B(D S -> tv) = (6.8 ± 1.0)% , 

where the uncertainties are statistical only. 

6 Conclusions 

Leptonic decays of the D s meson have been studied in a sample of four million hadronic 
Z decays. The D s — ► tv analysis gives 

B(D S -> tv) = (5.86 ± 1.18 ± 2.16)% , (electron channel) 
= (5.78 ± 0.85 ± 1.81)% , (muon channel) 

where the measured signal includes some D s — > fiv and D + — > Iv decays. In the Standard 
Model, these branching fractions correspond to the following values of the decay constant 
/ds : 

/d. = (275 ± 28 ± 51) MeV , (electron channel) 

= (273 ± 20 ± 43) MeV . (muon channel) 

These results are combined, taking into account the common systematic errors |24}| , 
to obtain 

B(D s ^tv) = (5.79 ± 0.77 ± 1.84)% 
/ Ds = (273 ± 18 ± 43) MeV . 
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The combination has x 2 = 0.0024 for 1 degree of freedom (CL = 96%). This measurement 
may be expressed as 

f(c -> D S )5(D S -> tv) = (6.77 ± 0.90 ± 1.49) x 1(T 3 , 

where f(c — > D s ) is the D s production rate per hemisphere in Z — > cc events; the relative 
uncertainty on this product is smaller than that on £>(D S — > tv) due to the reduced 
dependence on the assumed D s production rate. 

The analysis optimized for D s — > fiv is also sensitive to D s — > tv decays; a measurement 
of the combined signal yields 

B(D S -> fiv) = (0.68 ± 0.11 ± 0.18)% 

when the ratios of the signal components are fixed to their Standard Model values. This 
measurement corresponds to 

/ Ds = (291 ± 25 ± 38) MeV . 

Finally, the three /□„ measurements (two from tv and one from fiv) are combined. 
The correlation coefficient k = +0.43 ± 0.34 of the statistical errors in the D s — > fiv 
analysis and the muon channel of the D s — ► tv analysis is evaluated by dividing the data 
into ten subsamples. The correlations between the systematic errors are also taken into 
account. The combined decay constant is 

/ Db = (285 ± 19 ± 40) MeV 



with x — 0.51 for 2 degrees of freedom (CL = 78 ( 



This result is consistent with recent measurements ||T6| , |25| , |26|| and with the lattice 
QCD prediction, / Ds = 255 ± 30 MeV Q. 
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